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Abstract 

The simplest unified extension of the Minimal Supersymmetric Standard Model 
with bi-linear R-parity violation naturally predicts a hierarchical neutrino mass 
spectrum, suitable to explain atmospheric and solar neutrino fluxes. We study 
whether the individual violation of the lepton numbers L e fl<T in the charged sec- 
tor can lead to measurable rates for BR(/i — > 67) and BR(t — > ^7). We find 
that some of the R-parity violating terms that are compatible with the observed 
atmospheric neutrino oscillations could lead to rates for [i — > ej measurable in pro- 
jected experiments. However, the Amf 2 obtained for those parameters is too high 
to be compatible with the solar neutrino data, excluding therefore the possibility of 
having measurable rates for \jl — > ej in the model. 
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1 Introduction 



In the Standard Model (SM), lepton number is exactly preserved in contradiction with 
the observed neutrino oscillations |[[ |2|. Several extension of the SM include patterns of 
neutrino masses and mixings which can provide a satisfactory explanation for these flavor 
oscillations. The consequences of the individual violation of the lepton numbers L ejjUjT for 
charged lepton will be manifest in processes such as fi — > ey, \x — > 3e, fi — e conversion 
in heavy nuclei, r — ► /ry and Kl — ► /xe ||. The experimental upper bound for these 
processes is quite restrictive, which imposes a significant constraint for the explanation of 
flavor in models beyond the SM. However, the mechanisms used to explain the origin of 
the tiny neutrino masses required to explain solar and atmospheric neutrino oscillations, 
typically imply that these processes may occur at small rates, motivating an increasing 
experimental interest in exploring further charged lepton flavor violating processes. 

The rates for charged lepton flavor violation (LFV) are extremely small in the SM 
with right-handed neutrinos (oc AmJ/M^, Q). In R-parity conserving supersymmetric 
(SUSY) models, like the Minimal Supersymmetric Standard Model (MSSM), the presence 
of LFV processes is associated with vertices involving leptons and their superpartners || . 
These processes are sensitive to the scalar mass matrices structure, a non-diagonality of 
the latter in a basis in which fermions are diagonal, leads to a hard violation of flavor. The 
structure of the scalar mass matrices is very sensitive to the SUSY-breaking, in particular 
in models where SUSY is softly broken, LFV imposes a severe constraint in the flavor 
dependence of the soft-terms as they are generated in GUT's and string inspired models 

The inclusion of a "see-saw" mechanism in the MSSM provides a very attractive 
scenario to understand neutrino oscillations with very small neutrino masses, and at the 
same time gives rates for LFV processes accessible in projected experiments J?], [J. The 
waiving of the R-parity symmetry in the MSSM provides an alternative scenario to explain 
the generation of small neutrino masses. In this case the R-parity violating operators can 
be constrained by rare processes [§, [L(], [11], • 



The simplest extension of the MSSM with bilinear R-parity violation (BRpV) JO 
(allowing B-conserving but L-violating interactions) can explain neutrino masses and 
mixings which can account for the observed neutrino oscillations [|14j|. The BRpV model 



has been extensively discussed in the literature fL6 |. It is motivated by the fact that it 
provides an effective truncation of models where R-parity breaks spontaneously by singlet 
sneutrino vev's around the weak scale [[17]] . Moreover, they allow for the radiative breaking 



of R-parity, opening also new ways to unify Gauge and Yukawa couplings |18| and with 



a potentially slightly lower prediction for a s ]B| . For recent papers on phenomenological 



implications of these models see Ref. [20, 21]. As the parameters involved in the R-parity 
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violating operator are constrained in order to predict neutrino masses in the sub-eV range, 
we address in this paper the question of whether this operator will induce rates for charged 
LFV processes of experimental interest. Some of them occur at tree-level such as double 
(3 decay (12], |22| and /x — e conversion in nuclei |p3| . One loop LFV decays as lj — > 
become interesting on this framework due to the experimental interest in improving the 



current limits |24j] 



BR(fi -> erf) < 1.2 x KT 11 
BR(t^hj) < 1.1 x 1(T 6 

BR(t -> ej) < 2.7 x 1(T 6 . (1) 

As we will show, the predictions for the last two processes are much lower than the above 
limits and will not constrain the BRpV model. For /x — > cy the predictions are compatible 
with the current limit but could begin to constrain the model for the bounds that will 



be reached in current |25[] or planned experiments [p^] , if only the atmospheric neutrino 
data were taken in account. However the requirement that the one-loop induced Amj 2 is 
in agreement with the solar neutrino data will imply that the predicted rates for /x — > e^j 
will not be visible, even in those new experiments. 

This paper is organized as follows. In sections [| and f| we describe the model, the 
scalar potential and the fermion mass matrices, respectively. In section |5| we derive the 
expressions for the LFV processes. The results are presented in section ^| and in section [7] 
we give our conclusions. The more technical questions regarding the mass matrices, 
couplings and the explicit formulas for the amplitudes are given in the appendices. 



2 The Superpotential and the Soft Breaking Terms 

Using the conventions of Refs. [^T|, we introduce the model by specifying the super- 
potential, which includes BRpV |16] in three generations. It is given by 

W = e ab [h^QtuM + h%Q h j5jR a d + h%l%E a d - ^H a d H b u + e^flj] (2) 

where the couplings hu, ho and He are 3x3 Yukawa matrices and \x and e« are parameters 
with units of mass. The second bilinear term in Eq. (||D violates lepton number and 
therefore also breaks R-parity. The inclusion of the R-parity violating term, though small, 
can modify the predictions of the MSSM. The most salient features are that neutrinos 
become massive and the lightest neutralino is no longer a dark matter candidate because it 
is allowed to decay. Furthermore, we can observe that this model implies the mixing of the 
leptons with the usual charginos and neutralinos of the MSSM. Lepton Yukawa couplings 
can be written as diagonal matrices without any loss of generality since it is possible 
to rotate the superfields L\ in the superpotential, Eq. (0), such that Yukawa matrix h E 
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becomes diagonal. Conversely, in BRpV models it is possible to apply a similar rotation 



to reduce the number e parameters and provide a non-trivial structure to He |28 



Supersymmetry breaking is parameterized with a set of soft supersymmetry breaking 
terms. In the MSSM these are given by 

Csoft = -V™ f f M + [±M S A A + \MXX + \M'X'X' + h.c] (3) 

where 

Ko/f M = M Q 2 QTQ a j + MflJiU* + MfDiD* + MfL^L" + MfRiR* 

+mj I H a d *H a d +m 2 Hu H:*H a u (4) 



+£ab 



AgQfUjHl + A%Q b i D j H^ + A^L^H* - B^H a d H b u 



In addition to the MSSM soft SUSY breaking terms in V S ^ SM the BRpV model contains 
the following extra term 

where the Bi have units of mass. 



VfoT = -B^eJJlEl , (5) 



The electroweak symmetry is broken when the two Higgs doublets H d and H u , and 
the neutral component of the slepton doublets Lj acquire vacuum expectation values. We 
introduce the notation: 

*-(£)• *-©• M?)' 

where we shift the neutral fields with non-zero vev's as 

H d = ^[a d + v d + i<p d ], fl2 = -L[ ff 2 + « B + .VS], L° = -^[g* + + (7) 



Note that the W boson acquires a mass = \g 2 v 2 , where v 2 = v d + v% + vf + v\ + v\ r - 
(246 GeV) 2 . 



In addition to the above MSSM parameters, our model contains nine new parame- 
ters, 6i, Vi and Bi. The minimization of the scalar potential allows to relate some of these 
free parameters. The values of e^, V{ are directly related with the neutrino masses and 
mixings as we will discuss below. 



3 The Scalar Potential 



The electroweak symmetry is broken when the neutral Higgses and the neutral slepton 
fields acquire non-zero vev's. These are calculated via the minimization of the effective 
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potential or, in the diagrammatic method, via the tadpole equations. The full scalar 
potential at tree level is 

2 



v° 

v total 



E 



dW 



dzi 



+ V D + V. 



MSSM 
soft 



BRpV 
soft 



where z% is any one of the scalar components of the superfields in the superpotential in 
Eq. (H), Vd are the D-terms, and V^ f^ v is given in Eq. (H). 



The tree level scalar potential contains the following linear terms 



V, 



o 

linear 



(9) 



where the different t° are the tadpoles at tree level, their explicit expressions can be found 
in Ref. The five tree level tadpoles t° a are equal to zero at the minimum of the tree 
level potential, therefore we can use them to express the parameters: 



//, B, Bi ,B 2 , B : 



in terms of: 



V u ,V d , 6i, Vi, M 

We have two possible solutions for \i: 



3- 



2 2 2 

Lip m H u i m H d - 



(10) 

(11) 



-b ± Vb 2 - Aac 



2a 



'12) 



where a = v% - v 2 d , b = 2v d £? =1 e^i and 



™>H U ~ V l m H d 



- E 



3 3 



-^-^1^ + ^) (13) 



where we have defined D 



i=i j=i 

2 „,2N 



;(r + g' 2 m + v\ + u§ + v\ - v 2 u ). 



As one can easily verify, the above relations lead to the MSSM relation for \j? in the 
limit of vanishing and t>j. The uncertainty of the sign in the MSSM is translated here 
into two possible values for /x-term. However for the values of 6j and v i relevant to our 
work both solutions are close in module and of opposite sign. The values for B and Bi 's 
can be expressed in terms of \i as: 

' (14) 



B = — 

v.. 



v d i 2 , 2 

— {m H + fi 

H d 



B; 



1 



1 



8=1 



0^0 



(Dv i + -j:v j (Ml ij + Ml ji )) 



(15) 



The equivalent equations for the MSSM equations are obtained by setting ej, Vi equal 
to zero. 
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4 Fermion Masses with BRpV 



As we discussed in the previous section the presence of BRpV terms in the superpotential, 
Eq. (0), induces non-zero vev's for the sneutrinos and enables the neutrinos to have a 
mass, with a value related to the size of e«, t>j and the SUSY parameters involved on 
the electroweak symmetry breaking. Furthermore, the non-conservation of the R-parity 
allows the SUSY partners to mix with the SM particles. In this section we describe with 
detail the resulting neutralino-neutrino and chargino-charged-lepton mass matrices, since 
they are the most directly related to our problem. The complete set of mass matrices for 
the BRpV Model can be found in Ref. 0. 



4.1 Neutralino— Neutrino Mass Matrix 



The range of values of the e-parameters is indirectly associated to the size of the neutrino 
masses predicted by the model. To explore this relation we describe next the mass mixings 
among neutralinos and neutrinos. In the basis ip 0T = (— iX', —i\ 3 , H®, H®, u e , u^, u T ) the 
neutral fermion mass matrix is given by 



where, 



M N 



Mi 


- \g'v d 



m 



M 2 

\gvd 
-\gv u 



-\g'v d \g'v u 
\gvd -\gv u 



o 







(16) 



(17) 



is the standard MSSM neutralino mass matrix and 

-\g' v i \g v i o e i 

m = - \g'v 2 \gv 2 e 2 
- §5^3 Igva e 3 
characterizes the breaking of R-parity. 

The mass matrix M \ is diagonalized by 

M*M N M~ 1 = diag(m x o,m Uj ), (19) 

where (i — 1, • • • , 4) for the neutralinos, and (j = 1, • ■ ■ , 3) for the neutrinos. 

Since m u m x o the mass matrix is similar to the "see-saw" mass matrices 



5 



and takes approximately the form diag(.M x o, m e //), with 



/ A 2 A e A, A e A T \ 
A e A^ Al A M A r 
\A e A T A M A r A\) 



x rp M ig 2 + M 2 g' 2 
m eff = -mM m = ...... ? 

x 4 det(A^ x o) 

where the Aj parameters in Eq. fl20| ) are defined as: 

Aj = fIVi + V d €i 

One of the neutrino species acquire a tree level non-zero mass, given by: 

*»* ~ Tr ( m eff) 



(20) 



(21) 



, , M l9 2 + M 2 g' 2 lXl2 

= l r[m P jf) = . - , . . — — |A| , 



4 det(A4 x o) 



:22i 



where |A| 2 = X^ =1 A 2 . The two other neutrinos can get masses at one-loop as it is 
discussed in Ref. [H]. For our purposes it will be important to have an estimate of the 



values of AmJ 2 = m^ 2 — m^. We will use the results of Ref. [15] where it was found that, 



to a very good approximation, m ui = and 



3 hi m i 

m v2 = ttt^ m b sin 26 b — log 

lo7r z u z mf 

^ bi 



(ex A) : 
|A| 2 



(23) 



The explanation of the data on neutrino oscillations given in Ref. [I3| requires the 
neutrino masses to be be on the sub-eV range in order to fit the data on atmospheric 
neutrino oscillations. In our examples we take a m U3 = 0.1 eV, which leads to values 
of the | A] in the range of 0.1 — 1 GeV 2 , for the values of the SUSY parameters that we 
will consider. Considering that we take positive values for /i we should also take negative 
values for the product to avoid our analysis to be constrained to small values of tj. 
However, as we will see, for the values of the SUSY parameters that give the largest 
BR(fi — > cy), the values of the |ej| have to be below 0.1 GeV, if Am 2 2 < 10~ 4 (eV) 2 , as 
required by the solar neutrino data. 



4.2 Chargino— Charged Lepton Mass Matrix 

Due to the R-parity violating terms in the superpotential, Eq. (0), the charginos mix 
with the charged leptons, linking therefore the problem of the masses of the neutrinos 
with the problem of the charged lepton flavor violation. We describe in this subsection 
the chargino-lepton mass matrix to explain how the flavor mixing on the charged lepton 
sector arises. In a basis where 

*P +T = HA + ,^+,e+,/i+T+) and = (-zA~, Hj, el, m, r L ) (24) 
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the corresponding charged fermion mass terms in the Lagrangian are 

I, l+T T ,( Mq \ ftp 



where the chargino-charged lepton mass matrix Mc is given in the Appendix [A|. As in 
the MSSM, Mc is diagonalized by two rotation matrices, and we include the physical 
charged leptons and charginos into a set of five charged fermions defined as: 

Fr = U ij ^J ; F? = V l3 ^ (26) 

Such that 

U'McV ' = M CD (27) 
where Men is the diagonal charged fermion mass matrix. 

In the previous expressions the F i are two component spinors. We construct the 
four component Dirac spinors out of the two component spinors with the conventions^ 

*r - (J I (») 

The parameterization of the matrices V, and U given in Appendix |A], that was 
introduced in Ref. |22|, |29fl , provides a very accurate representation of the exact result. 



By comparing the numerical results with the analytical expressions shown in Appendix 
A, we found discrepancies of less than the 1%. To obtain this level of accuracy we had 
to introduce corrections in the definition of Vl, Vr, and Qr with respect to the formulas 
of Ref. P2"| . These arise mainly from including the sub-matrix E' in our derivation (see 
Appendix A). Although the size of the matrix elements of E' is smaller than the other 
components of Mc, it must be taken into account in order to match the results of the 
smaller elements of U and V found in the exact diagonalization. Our definition of Vl and 
Vr leads to the correct form for the lower right 3x3 sub-matrices of U and V. We will 
make use of it to explain the details of our results. The inclusion of the matrix E' in 
the determination of Qr allows to display the dependence of the matrix elements on the 
A-parameters, rather than a explicit dependence on the e's as quoted in Ref. pT 



Also, we must observe that the elements of Ql exceed the ones of Qr by several orders 
of magnitude. Therefore we can anticipate that the couplings containing the matrix V in 
Appendix 0, will be suppressed with respect to the ones containing the elements of U . 



*Here we depart from the conventions of Ref. |2j because we want the e , fi and r to be the 
particles and not the anti-particles. 
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5 lj — > /j7 Flavor Violating Processes and the fi Ano- 
malous Magnetic Moment 



5.1 Effective Lagrangian and Diagrams 

The effective operators that generate the decays lj — > l^'j and the lepton anomalous 
magnetic moment can be written as: 




Figure 1: Generic Feynman diagrams for Ri j. 

The one-loop contributions to Al,r in the model under consideration arise from the 
diagrams of Figs. |T]-^| 

A Lij = A s LlJ + A% + A% (30) 

Amj = A Rij + A%j + A Rij (31) 

The partial contributions on the above expression correspond to the addition of the sets 
of diagrams represented in each figure: 

A G _ A N»-W± , A C±~Z° (n^ 

A s _ A N a ~s ± , a c± ~ s ° 4- A c± - p ° 

aQ _ Adrf-u | Ad-vq , A wy-d . AU-d-y 

A L,Rij — A L,Rij ^L,Rij ^L,Rij ~t~ ^L,Rij 

The superscripts in each contribution on the right denote the fermion and boson internal 
lines of the corresponding diagram. For the quarks- squarks diagrams we include the 
symbol 7 to indicate whether the photon is attached to the fermion or the boson line. 
We follow the notation of \T%\ indicating by S the eigenstates of the charged scalar 



mass matrix, by S° and P° the eigenstates for the sneutrino-Higgs scalar mass matrices, 
CP-even and CP-odd, respectively. 
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The contributions to Al Ri - arise from the diagrams in Fig. [I]. The index A = 1, ... ,5 
corresponds to the eigenstates of the chargino-lepton mass matrix, while the indices 
i,j = 1,2,3 correspond to the lepton generation indices in the limit of the MSSM with 
R-parity conservation. These diagrams will become the SM contribution to £ e //, Eq. (|29|) , 
in the limit q = 0. In the case of the SM it provides the main contribution to a^, no 
contribution for charged LFV processes when neutrinos are considered massless and a very 
suppressed contribution for the values of m Ui compatible with the experimental limits . 

The contributions to A£% arise from the three diagrams of Fig. ||], where the index 
X refers to the eigenstates of scalar mass matrices. S are the eigenstates of the 8x8 
charged Higgs-slepton mass matrix and S°, P° represent the eigenstates of the 5x5 neutral 
Higgs-sneutrino scalar and pseudoscalar mass matrices, respectively. In the limit where 
R-parity is conserved these three diagrams will be combined in the two supersymmetric 
diagrams contributing to the Al/r in the MSSM. In this this limit, these diagrams are 
flavor conserving when the soft-terms are universal as given by minimal Supergravity 
version of the MSSM. 



y 



S° (P°) 



X,\ 



Xi 



5Ca 



Xa 



xj 



y 



Figure 2: Generic Feynman diagrams for A£ l i R . 

A®/ R arise from the four diagrams of Fig. [| where the indices X = 1 , . . . , 6 refer 
to the eigenstates of 6 x 6 squark mass matrices and indices a = 1,2,3 are the quark 
generation indices. These diagrams are not present when R-parity is conserved. 



5.2 BR(lj — > ^7) for Flavor Violating Processes 

The branching ratio for the rare lepton decays lj —>■ 1^ is given in the literature 0] and 
we do not repeat the derivation here. The result is 

BR(IJ - Kl) = ^ (\A Lij \ 2 + \A RlJ \ 2 ) (35) 
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where the amplitudes A Li j and A Ri j were defined in Eq. (|30D . The complete expressions 
for the amplitudes corresponding to these processes in the BRpV model are given in the 
Appendix y. In their derivation we have neglected the mass of the outgoing fermion. 



5.3 The Muon Anomalous Magnetic Moment 



The expression for the muon anomalous magnetic moment can be obtained from the 



Lagrangian given in Eq. (|29|). One obtains |30| 



(36) 



The amplitudes A^j R can also be obtained from the formulas of Appendix C by including 
the effect of in both external lines of the diagrams. To do that we just have to include 
a factor of 2 in the part of the amplitude containing the function fp, P = N,C, W, Z: 



L/R 



Af /R (fp(x) 



2/p(x)) 



(37) 



y 

n x (3x)> n x (3 x ) 



d x (u x ) 



d a (u a ) 



X, u a (d a ) > u a (d a ) % 



Figure 3: Generic Feynman diagrams for AfJ, R . 



6 Results 



6.1 The Parameter Space 



The BRpV model that we consider adds more free parameters to the ones already present 
on the MSSM. However, if we consider the phenomenological constraints imposed on 
the MSSM by the limits on the mass of the lightest neutral CP-even Higgs m^, by 
the BR(b — > 57) and by the value of the a M , as well as those derived from neutrino 
physics on the BRpV parameters, we can narrow the space of parameters such that 
generic predictions for BR(fi erj) and BR(t — > ^7) can be made. 
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We assume the parameter space of the MSSM with universal soft-terms and GUT- 
unification, 

«g, M GUT , mo, M1/2, tan/5, //, 5, ft^, A (38) 
with the addition of the BRpV parameters, 

e i7 Bi, Vi, i = 1,2,3. (39) 

A is defined such that Aj(GUT) = A - hj, I = U, D, E. The quantities otc = Qg/^ (#g 
being the GUT gauge coupling constant) and Mqut are evaluated consistently with the 
experimental values of a em , a s , and sin 2 ^ at rriz- We integrate numerically the RGE's 
for the BRpV model, at two loops in the gauge and Yukawa couplings and a one loop in 
the soft terms, from Mqut down to a common supersymmetric threshold Ms ~ ^ m t 1 m t 2 ■ 
From this energy to mz, the RGE's of the SM are used. 

As we explained before, the minimum conditions of the effective scalar potential 
allows us to express the values of /i, B, B\ B 2 , B 3 in terms of tan/3, 6j, Vi. These are 
evaluated at the scale M$. The value of fi obtained at this scale is similar to the one 
obtained by minimizing the effective potential with the complete 1-loop MSSM contribu- 
tions [[31]]. The 1-loop contributions arising from Rp- violating terms for these parameters 
are comparatively much smaller. 

We fix the elements of the quark Yukawa matrices at the GUT scale, consistently 
with the experimental values of the quark masses and the absolute values of the CKM 
matrix elements. In the case of the charged leptons we have to make sure that the 
three lightest eigenvalues of the chargino-charged lepton matrix are consistent with the 
experimental values of the charged lepton masses. 

The values of mo and are chosen in region of the parameter space favored by 



the considerations presented in Ref. [30|, so that we can compare our results with typical 
predictions for the BR(fi — > ej) in the MSSM with a "see-saw" mechanism, as discussed 
in Ref. M. Obviously, since our model breaks R-parity, the LSP is not a dark matter 



candidate and therefore the cosmological preferred areas of Ref. p0[ do not apply to our 
study. However, the restriction of considering points in the m — m 1( / 2 plane such that 
irih > 113 GeV, is the most restrictive. The SUSY contribution to a M favors the sign 
of fi to be positive for the choice of SUSY parameters given below. We found that the 
upper bound of Eq. fl48|) (see below) on 5a^ is less restrictive than the one imposed by 
m h > 113 GeV. We analyze three sets of SUSY parameters, 



a) tan/3 = 10, m 1 / 2 = 400 GeV, m 

b) tan/3 = 30, m 1/2 = 400 GeV, m 

c) tan/3 = 30, m 1( / 2 = 600 GeV, m 



200 GeV, A = 0,m v = 0.1 eV. 
300 GeV, A = 0, m v = 0.1 eV. 
300 GeV, A = 0, m v = 0.1 eV. 



11 



The six free BRpV parameters £», V{ reduce to three if we take into account the 
constraints imposed by the predictions for neutrino oscillations in this model, as given 
in Ref. [TJ[]. By setting the atmospheric neutrino anomaly scale to the magnitude of the 



tree level non-zero value of one of the neutrinos, Eq. (|22]) , we fix the value X^Af for 
each SUSY point, where the Aj were defined in Eq. ([H]). We then follow the discussion 
of Ref. ||14|| , where it was shown that the conditions A3 ~ A2 — 5 x Ai satisfy both the 



atmospheric neutrino anomaly mixings and the CHOOZ result 0. We then obtain a 
linear relationship between each couple e», Wj. 

Therefore we study the dependence of the process under consideration on the values 
of £1, €2, €3 for a neutrino mass of m U3 = 0.1 eV, on the upper limit of the allowed range 
for the atmospheric neutrino anomaly. For comparison purposes we also present some 
results for m„ 3 = 0.05 eV on the middle of that range, anf for m Vz = leV. 

We will assign random values to e\ and e 2 in the range: 

- 2 x 10~ 3 GeV > ex, e 2 > -60 GeV. (40) 

However the requirement that m U2 < 0.01 eV will exclude values of |e$| > 0.1 GeV. This 
will be explicitly shown in our results. The value of €3 is kept fixed since our results are 
not altered when it varies on the above range. 



6.2 The Branching Ratio for lj — > l { 7 

We perform a full numerical analysis with the exact diagonalization of matrices involved 
in the computation of branching ratios. The main contribution for BR(fi — > ej), Eq. flSS]) , 
comes from the amplitudes Ar. The partial contributions from the various diagrams listed 
in Eq. fl32"H34]) are displayed in Fig. [|, for the set of parameters b). We have found that 
they are all independent of €3 and that they display a linear behavior as a function of 
the product e\ ■ €2, when this product is larger than 0.1 GeV except for the cancellation 
observed in A^ = Ar ~ s ° + Ar ~ p ° . The values of the amplitudes arising from the 
diagrams of Fig. [1] depend on the e's through the A's which are kept fixed, and therefore 
remain constant. We also found the contributions of the diagrams of Fig. |3] to be of the 
same order of magnitude. As we can see, the sum of all of them, A^, is almost a linear 
function of e\ ■ e 2 . The amplitudes arising from the diagrams of Fig. |2] are the dominant 
ones. The one mediated by the neutralino (^4^) is smaller than the dominant chargino 
exchange (A^), except in the range where the cancellation takes place. 

We can also observe in Fig. |] that the cancellation which appear in the A^ depends 
on the value of e±. The values that we show correspond to two different choices of e\ (note 
that if we allow e\ to change randomly, as we have done with the other amplitudes, the 
values for A^ would be a distribution of dots). The behavior of A^ cannot be attributed 
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Figure 4: Partial Amplitudes Ar corresponding to the choice of parameters b) of section 
6.1 for ei = —5 GeV (solid) and e\ = —0.1 GeV (dash). The values of change sign on 
each branch of the curve, the left one corresponds to positive values. 

to an accidental cancellation between the scalar and pseudoscalar parts as one may naively 
expect, on the contrary both parts add constructively and almost vanish simultaneously. 
The behavior of that amplitude can be explained when we identify which are the particles 
running in the loops of Fig. |2]that are responsible for the main contributions: X — 1, A = 
1 and X = 4, A = 1, 2. Then we can obtain an accurate approximation for the amplitude 
by using the formulas given in the Appendices |A], || and |C| Let's consider v4^ ±_s ° since 
the contribution of the corresponding pseudoscalar exchange is almost identical. We get 
from Eq. (^) for the dominant contributions, 



A' 



1 



32tt 2 



I II go "'/I Si A=l A 1 



(41) 



Using the definitions of Appendix || we find, 



V^iiVlTi « U 32 U U , Vg^VfiU « * U U V A1 U A2 (42) 

where is the Yukawa coupling of the muon. We can then write, 

^m4~ S ° « *i U 32 U 14 + F 2 U 34 (43) 
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where 



gK 1 



F i = i^h;^r h c{xu)— (44) 



nh 1 2 

E^(^)— (4 5 ) 



The quantities iq and _F 2 are independent of the BRpV parameters and can be eval- 
uated given the SUSY parameters. The dependence of the amplitude Ar\a~ s ° on the ej 
comes from the matrix elements U32 , Uu and t/34. Using the Appendix [A] we can find ap- 
proximate expressions for these matrix elements that display explicitly this dependence. 
We get, 

U M ~ r + 2, ^32 ~ — , 2/i4 ~ -Ui 2 —. 46 

A* "^d/ 1 A* A* 

Hence we can find the value of e 2 at which Ar ~ 0, 

62 ~ (47) 

The position of the cancellation changes, with the value of the SUSY parameters and also 
with value of A 2 , as we can see from Eq. fl4"T|) . This explains the qualitative changes we 
find in Figs. 



Some of the the amplitudes contributing to the BR(/i — > ej) presented above have 
been previously discussed in Refs. || pTTj . We agree with Ref. [JO]] in that the main 
contribution arises from except for the values of parameters affected by the cancellation 
mentioned above. However we find smaller values for A® than the ones quoted in || . 

The contribution of Ai to the branching ratio is negligible compared with Ar, due 
to the fact that the matrix U is replaced by V with suppressed mixings. This holds even 
for the element V34. As we can see in Appendix |A] this element is determined by Vr which 
is obtained in a similar way as Vj, for the matrix U. However, we observe that the main 
contribution to V34 is suppressed by a factor m e /m^ with respect to the corresponding 
one in U34. 

Fig. [| shows the impact of the cancellation in A^ on the predictions of BR(/i — > ey) 
for the choice of parameters b). As we can deduce from Eq. (|47[) the value of the e 2 at 
which the cancellation in A^ takes place depend on the values of the SUSY parameters. 
This determines the shape of the curves of constant BR in Figs. ^ |7| and |8]. Since the main 
contribution comes from the chargino mediated diagram of Fig. 0, we can expect that the 
BR increases with tan/3 and decreases as m^j grows. The increase of mo produces the 
same qualitative effect as the increase of mi/ 2 , however it has a lower impact on BR than 
the changes in m^j. 

Our results can be compared with the predictions of a model based in the MSSM 
with a "seesaw" mechanism presented in Ref. ||, where the results for BR{^i — > ej) are 
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of order 10~ 13 for case a) and between 1CT 12 and 10~ 13 for b) and c). If we observe our 
predictions for case b) on the left graph of Fig. || we can see that the model predicts 
ratios of 1CT 11 — 1CT 13 for values of |ei| and | e 2 1 ranging from 1 to 10 GeV (independently 
of the value of €3). Values in the range of 0.1 to 1 GeV would lead to rates of order 
10 -14 — 10 -16 , still interesting for the next generation experiments |25l [26|]. Similarly, a 



window of 0.1 < — e 1; — e 2 < 1 GeV, is crossed only by the 10~ 16 line in case c) and by 
lines below this value for case a). Such values of \e{\ are however excluded if one takes 
in account the constraint coming from the solar neutrinos mass scale. This is shown in 
Figs. |6| and |7], where the dashed line gives the upper limit on the |ej| as obtained from 
Eq. ( p3[ ) for the requirement that m U2 < 0.01 eV . 



In these figures we can appreciate that the parameters that enhance the ratios also 
make m U2 larger. From Eq . (^3|) we can observe that m V2 increases with tan (3 (through 
the dependence on hb) and decreases as the /x-term increases (i.e with mi/ 2 and m ). 
Furthermore, since we have choosen the A-parameters to be related, m„ 2 is proportional 
to a combination of ef with no accidental cancellation amoung them. Therefore we can 
not find suitable values for tan/3, m\/2 and mo, such that we can find an overlaping of 
areas with m U2 < 0.01 eV and BR(/i — > ej) > 10~ 16 . 

In Fig. H we consider m U3 = 0.05 eV and m V3 = 1 eV both for case b). This decreases 
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-Ej (GeV) -Ej (GeV) 

Figure 6: Contour plot for BR(/i — > e7) in 61-62 plane for case b (left) and c (right). The 
dash lines correspond to m v% = 0.01 eV. 



(increases) the values of the A's by about a factor of \/2 (vlO), respectively. By looking 
at the parameterization of the matrix U in Appendix A we can infer that these changes 
in the A's have not a decisive impact on the /1 — > rates. The reason for this is that the 
dominant contributions to A R34 ' are determined by the matrix Ql and its elements 
that depend explicitly on the e's are much larger than the ones containing A's (at least 
for values of e's leading to relevant ratios). However the position of the cancellation on 
^i?34~ S ° depends on A 2 as we can see in Eq. (|47D and it therefore determines the changes 
in the figures. 

The changes on the A's have only a direct impact on the smaller contributions, such 
as the ones arising from the diagrams of Fig. [I] and on the Al, which size is controlled 
by the elements of Qr, which, as we have said, are several orders of magnitude below the 
main contribution coming from Ql- 

The predictions for BR(t — > /ry) that we obtain with this model are of the same 
order as those for BR(fi — > e^), the reason being that we have assumed the A's to be of 
the same order of magnitude, as it is required to explain neutrino oscillations. The results 
in this case are independent of e\. If we consider values for e 2 and €3 in the same range 
as in Fig. ^| we obtain similar curves. This result contrasts with the LFV results on the 
framework of the R-parity conserving MSSM, where BR(fi — > ej) is typically suppressed 
by several orders of magnitude with respect to BR(t — > yrf). In this case the hierarchy 
of Yukawas couplings makes a distinction between the two processes. 

We conclude therefore that the values of the parameters of the BRpV model that 
successfully explain the solar and atmospheric neutrino data [14]], predict rates for // — > e^y, 



fij and t — > ej, that are well below the current limits and well as those of planned 
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Figure 7: Same as Fig. for case a. 



experiments. 



6.3 The Muon Anomalous Magnetic Moment 



The difference on the value of the muon anomalous magnetic moment found in the BNL 



E821 measurement [S3| with respect to the SM prediction, which originally was considered 



to be 2.6 a is now reduced to 1.6 a after a theoretical error has been corrected [ 34[ . When 
the 2 a range is considered, the allowed values for contributions beyond the SM become, 



- 6 x 1(T 10 < 5a 



exp 



a s u M < 58 x 1(T 10 , 



(4£ 



Several studies 30 indicate that the MSSM extension of the SM can account for this 



discrepancy. When R-parity is broken the SUSY particles are allowed to enter in the SM 
diagrams (Fig. [I]) and conversely the SM particles run in the SUSY loops (Figs. |2] and |3[). 

The contribution due to the R-parity violating operators to 5a^ is obtained by 
subtracting from the amplitudes arising from Fig. [TJ (da^(RpV), 5a^ (RpV)) and Fig. ^] 
(5a% (RpV), 5a*° (RpV)) the contribution of the SM and the MSSM, respectively (which 
are obtained in the limit of vanishing e^'s). The contribution from Fig. |3] (a^) is not 
present in the MSSM. All these contributions are found to be small when the R-parity 
violating terms are associated with neutrino masses of experimental interest. 
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Figure 8: Same as Fig. || case b for m U3 = 1 eV (left) and m V3 = 0.05 eV (right). 









a q 


oaf{RpV) 


5af(R P V) 


SaZ(RpV) 


Sa^(RpV) 


a 


9.6 


0.45 


-1.5- 10~ 4 


7.5 • 10" 2 


6.1 • 10" 3 


-0.25 


0.51 


b 


25.8 


-1.1 


-1.5 ■ 10" 3 


0.35 


2.5 • 10" 2 


-0.28 


0.52 


c 


13.8 


0.19 


-4.1 • 10" 4 


0.11 


1.9 • 10" 2 


-0.14 


0.28 



Table 1: Contributions to a M from the graphs in Fig. []V|3| (in units of 10~ 10 ). Cases a-c, 
refer to the choice of parameters given in section |6J]. See section [5] for details on the 
notation. 



In Table [l] we show the different contributions to 5a^ for the selected values of 
m 0) m i/2 and tan (3 discussed in the section p]T|. The main contributions to 8a ^ arise basi- 
cally from the MSSM components of the diagrams in Fig. (a^ + , &?f). The contribution 
from BRpV operators just adds a small percentage to the total values arising from physics 
beyond the SM. The values that we show correspond to the maximum value obtained in 
the conditions for the e-parameters described in section |6.1| , when we allow |ei|, \t2\ to 
range from to 60 GeV (the result is almost independent of the value of e 3 ). 



7 Conclusions 

We studied the LFV in one-loop induced rare processes lj — > 1^, % ^ j in SUSY models 
with bi-linear R-parity violation. In this context, the R-parity violating interactions can 
explain the neutrino masses and mixings without adding new fields to the particle content 
of the MSSM which represent an appealing alternative to the "see-saw" mechanism. In 
this work we addressed the question of whether neutrino masses in the sub-eV range can 
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be compatible with rates for charged LFV processes of experimental interest. 

We have performed an exhaustive study of the interactions and of the SUSY param- 
eters involved in the processes lj — > lij, i 7^ j and contributing to the muon anomalous 
magnetic moment. For the case of rare decays we find the diagram mediated by the 
Higgs-sneutrino scalar s to be the dominant. Contributions arising from diagrams in- 
cluding lepton-squark and lepton-gauge bosons vertices, possible in this model, are very 
suppressed in the range of neutrino masses considered in this work. Regarding a^, the 
additional contributions introduced by the R-parity violating interactions modify in a 
small percentage the value obtained in the MSSM limit. 



As in a previous analysis [|nj we find the BR(/i — > cy) to be very sensitive to the 



product ei -t2- However the presence of a cancellation in the main amplitude contributing 
to this process (which we have analyzed in detail through an accurate parameterization of 
the matrices U, V), makes our contour plots sensitive to the values of e\ and 62 in most of 
our examples. The rate increases with tan /3 as it is the case in the MSSM with "see-saw" 
mechanism. However as the one-loop induced neutrino masses will also grow with tan f3, 
the requirement that Amf 2 is compatible with the solar neutrino data excludes the region 
in parameter space where the BR(/i — > e^y) could be of experimental interest. On the 
other hand the rates for r — > /i7 found in our study are of the same order as the ones for 
— * cy, therefore also out of the experimental range. 

Unlike the situation in these models, the rates for r — > /i7 found in our study are of 
the same order as the ones for /i — > ej, therefore out of the experimental range. 

To conclude, we must say, that the obtained results for the /i — > e'y show us that if the 
BRpV model is the explanation for both the solar and atmospheric neutrino ocillations, 
the predicted LFV will not be testable at PSI |25[ or at PRISM |p6|| . The correlations of 



the BRpV parameters with the neutralino decays, as proposed in Ref. [pq] , will remain 
the main test of the model. 
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A Chargino— Charged Lepton Mass Matrix 



The chargino-charged lepton mass matrix, in the basis of of Eq. (|24|) , takes the form 



M, 



c 



M x ± E' 
E M E 



(49) 



where Me = -j=Od\iE is the charged leptons mass matrix and M x ± is the usual MSSM 
chargino mass matrix, 



The sub-matrix E is 



E 



M 2 -^gv u 

^gv d /i 



T2 gV2 ~ 62 
j= 2 gv 3 -e 3 



(50) 



(51) 



and E' can be written as E' = —v. He, where v is defined as: 





V\ t>2 t>3 
V2 V? 



(52) 



As the R-parity breaking parameters are small compared with the SUSY scale, 
it is possible to have an approximate diagonalization of Mq- This will be very useful 
in understanding the numerical results as we can have approximate analytical formulas. 
This approximate diagonalization is obtained by using the following parameterization, 
introduced in Refs. [H, HJ, for ^*and 



U* 



ut 
o V L 



-n T . 



Qr 1 — ^QrQ r 



R 



\n L o} L 



u R o 



V, 



R 



where Ulr are the MSSM rotation matrices, 



U* L M X± U R 



M d '+ 9 



(53) 



(54) 



(55) 



The matrices £Il,r and Vl,r are to be determined from the unitarity of U and V, and 
from the defining condition 



U*M C V ] 



M diag 
X 

o 



o 

M diag 



(56) 



In the literature [22, 29 1, the matrices £Il,r and Vl,r were obtained in the approxi- 
mation E 1 = 0. However we discovered that this approximation was not good enough to 
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explain our numerical results. So we re-derived the expressions for these matrices without 
neglecting E' . We get the following expressions for Ql,r, 

i(ei det(M x± ) + |sVAi 

&l = EM;l = - ^ , ^A 2 -i(e 2 det(M y± ) + ±<? 2 ^A 2 ) (57) 



det(M 



^A 2 - J(e 2 det(M x± ) + |^ 2 ^A 2 ) 
^A 3 -i(e 3 det(M x± ) + ±<r^ u A 3 ) 



n 



R 



(E« + M\EM~l) (M~l) T 


= h E \-v 


f + 








iAi 


1 

det(M x ±) 


-^(M E ) t2 A t 




2 Ai 




_ ^{M E )^A t 


~^ E ) 


3 Ai 



r-l\T 



-1\T 



(5? 



where summation over i = 1, 2, 3 is implied in each matrix element. The expression for 
Ql coincides with the one found in the literature but the expression for flu it is different. 
For Vl r we found that instead of the relation pi 



V L M E V ] R = M d E iag 



(59) 



they should satisfy, 



v L (m e - n L E' - \n L n T L M E ) v ] R 



M diag 



(60) 



For a general form of the matrix M E it will be difficult to have an explicit form for Vl,r- 
However for the case, that we consider, where the matrix M E is diagonal, we can obtain 
an analytical approximate expression for these matrices, 





1 


L.R 
Vl2 


L,R 
Vl3 


V L ,R ^ 


*L,R 
~ Vl2 


1 


L,R 
V23 




* L.R, 
~Vl3 


*L.R 
-V23 


1 



(61) 



where 



4 



€i€j_ m 2 + m) 

2 
3 



2 \i 2 mf — m 



+ 



9 v u v d 



m A 



9 v u v d 



2 det(M x 



2 det(M x 



2 2 

mf — rrij 



2 2 

— rrij 



9^Vu ™% + 

det (M x ±)/i 2 Vd mf — m 2 



AjAj 



(62) 



€i€j VTiiVTLj 



€jAj miTTlj 



fi 2 mf — m 2 fi 2 Vd mf — m 2 fi 2 Vd mf — m 2 



2 g 2 v u mimj 



det(M x ±)/x 2 Vd mf — m 2 



(63) 
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and rrii are the charged lepton physical masses. 

Putting everything together we can find analytical expressions for the matrix U that 
will be useful in explaining our results. We get, 

9 A 

<x 



u;+„2 - -n™-- (65) 

U* 1>2+i * U* Ll2 n Li2 ~-U* L12 6 -i (66) 

DJ.a+i - ^22 n M ^-^22 ^7 (67) 

/i 

UZhm ~ {Vl-I^V).. (68) 

For further reference we give the approximate expression for U^ A . We get 

U* e x e 2 £l A 2 g ^e 2 Ax 

34 " ^ 2 + /i 2 t; d det(M x± )/i 2 1 j 

- -^ + ^ (TO) 



where we have assumed that the parameters are in the ranges described in section 6.1. 



B The Couplings 

The relevant part of the Lagrangian, using four component spinor notation is, 
£ = [^(VmxPL + V^xPR)x° A Sx + ^(V^ x PL + V^ x P R )d A u* x 

+ xtivdxPL + v£ x P R )u A d x + xfnV^rPL + V^P^XaW- + h.c. 

+ xf (VETaxPl + v c R t AX p R ) XA s x + xT(V% p AX Pl + V^ AX P R ) XA Px 

+ X^r{V c L f A P L + V^P R ) XA Zl (71) 

The definition of these couplings is given in the following sections. These definitions 
extend those of Ref. [fill which conventions we follow. 



B.l Chargino-Neutralino- Charged Scalars 



■trans _ 
V Li A X - -PA 



gR s X2 (^N A2 V* 2 +N* M V^j +g' (^R&N^VZ +V2R s X5+a V i 
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VrTax = Vi gR S xi (A^N A2 U i2 - N A3 Ui^j + gR x % +a (-^U i2+a N A2 - N M+a U^\ 

+ (^f'i^l^i2+^X2+a^i2+Q^ll)+(^i2^V i 14+a — N A3 Ui 2+a ) Jljf R X 5-\ 



(72) 



(73) 



where the indices have the following ranges: A = 1, . . . , 7, % = 1, . . . , 5, a, ft = 1, . . . , 3 
and p A (rji) are the signs of the neutralinos (respectively charginos) as they are obtained 
from the numerical evaluation of the eigenvalues |14j . 

B.2 Chargino-Chargino-CP Even Neutral Scalars 



T/CCS 
V T A 



LiAX 



V2 



g (RxiUltaVn + R X2 U Al V i2 + R X2+a U A2+a V *l 



i2+f3 



Trees 
V RiAX 



V'CCS* 
LAiX 



(74) 
(75) 



B.3 Chargino-Chargino-CP Odd Neutral Scalars 



v LiAX 



IT)A 



V2 



g (R X iU A2 V^ + R X2 U A1 V* 2 + R X2 + a Ul 2+ a V ii 



+ (U A2 Rx2+a - R S X\U A2+ o) tifVi 



i2+/3 



TrCCp 

v RiAX 



V T CC F* 



LAiX 



(76) 
(77) 



B.4 Chargino-Neutralino-VK 



± 



V-cnW 
LiA 



-ty Pa g 



NI2UH + -j= (N* A3 U l2 + N* M+a U i2+c 



v; 



cnW 
RiA 



■j=N M V£ - N A2 V* 



(78) 
(79) 



B.5 Chargino-Chargino-Z c 



VccZ 
LiA 



ViVA 



cos9 w 



UnU* A1 + ( - - sin 2 9 W ) 5 iA 



JO) 



23 



V'eeZ 
RiA 



COS 9 



v^VAL + hfaVju-wPewSi. 



(81) 



B.6 Chargino- Quark Down-Squark Up 



V'cdu 
LiAX 



V'cdu 
RiAX 



„T/* D« pd* j T/* pi* !, Q P DM 

y v il rL Xa rL LAa "r v i2- n 'LAa' i ;7 rt X3+/3 



a/3 n» 



(82) 
(83) 



B.7 Chargino- Quark Up-Squark Down 



cud 
LiAX 



cud 
RiAX 



ViVA 



„ff* pa p«* tt* p«* !, Q /3 prf 
~9 u il n Xa n LAa "+* u i2 n LAa n D n X3+(3 



T/ P d (i^pi 
V i2 n Xa n U U RA/3 



U) 
15) 



C Amplitudes 

We collect here the various amplitudes corresponding to the diagrams of Figs. |]-|3|. In 
these amplitudes the mass of the outgoing fermion was neglected. We give only the 
amplitudes because the Ar can be obtained from these with the substitution rule 

A Rij = A Lij {L/R^R/L). (86) 



C.l Neutralinos— Charged Scalars 



aN°-S ± 
A Lij 



5 8 

EE 



i i 

32tt 2 ml ± 



A=l X=l 
( m x° V 

with %ax = ( ) and the functions fN>h N given by 



fN(x AX )V™Z x V™Xx + h N (x AX )^V^ x V%% x 



m ± 



17) 



fx(x) 



1 — Qx + 3x 2 + 2x 3 — 6x 2 In x 
6(1 -x) A 



h N (x) 



1 — x 2 + 2x In x 
(1-x) 3 



C.2 Charginos— CP Even Neutral Scalars 



A c 



-S° 



5 5 

EE 

A=l X=l 



32tt 2 



m 



fc(x AX )V£? AX V[^ x + h c (x AX ] 



ma 



'trees t^ccs* 
'LiAX v RjAX 



19) 



24 



with xax = y 

fc(x) 



m ± 



and the functions fc, he given by 
2 + 3x-Qx 2 + x 3 + Qxlnx 



6(1 -xf 



h c {x) 



-3 + 4x-x 2 -2\nx 
(1-xf 



C.3 Charginos-CP Odd Neutral Scalars 

1 1 



a=ix=i 327r 2 m p o _ 



fc(x AX )VlZxV^x + h c {x AX )^V^ AX V^X x 

mi, 



with xax 



m ± 



C.4 Quarks— Squarks 

~ 3 6 111 



. , v , 3 ; 32ttW 

A=l X =1 ux 



fc(x AX )ViSxV£f A *x + hc(x AX )^V^xV^ x 



with x AX 



u x 



~ 3 6 o i i 



A=1X=1 3 3271 " 2 m dx L 



~- — rri 

fc(x AX )V%i x V%f x + ^(^x)^^Ax^tx 



with aj^x 



4# 7 = E E 3(§) 1 : 



M*^)V£&V$k + h N (x AX )^V^xVgXx 



with = 



m dA 



^7/ 7 =E E 3 K) 



36 11 



A=l X=l 



3 y 327r 2 m 2 ~ 



}n{x ax ) Vmlx v SjAx + hx{x AX ) V^ AX V^ AX 



with x AX 



C.5 W-Neutralinos 



For W-Neutralinos the amplitude, in the unitary gauge (£ — > +oo) is, 



Lij 



E 



i 



f (m \l/cnWi/cnW* , /, /„ \ m X°A \rcnW '\rcnW* 

Jw{x A )V RiA V RjA + h w {x A )— — V mA V LjA 
J mi. J 
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with x A = ( ^7 ) and the functions fw, hw given by 



, , N 10 - 43a; + 78a; 2 - 49x 3 + 4x 4 + 18a; 3 In x . „. 

fw(x) = (9T) 

, . . — 4 + 15a; — 12a; 2 + x 3 + 6x 2 In a; . . 

h w (x) = -— -~ 3 (98) 



C.6 Z— Charginos 



For Z-Charginos the amplitude, in the unitary gauge (£ — > +oo) is, 



5 < j jyi ± 

C ±-Z° _ 1 1 f /„ \ V ccZ v ccZ* , i ( \ X A rrccZrrccZ* 

Jz{xa) v RiA v RjA + n z [x A ) V RiA V LjA 

mi 



Ll3 " h 32vr 2 m| 



(99) 



"m ± ■ 

with x A = ( ) and the functions fz, hz given by 



, . , 8 - 38x + 39x 2 - 14a; 3 + 5a; 4 - 18a; 2 lnx . . 

fz{x) = m 

h z (x) = (U ''" r (101) 
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